bet         = (1/(1.05))^(1./12) ;
N           = length(Month) ;
indvar      = linspace(1,N,N)-1 ;
options     = optimset('Display','off','TolX',1e-4,'TolFun',1e-4,'MaxFunEvals',1e10,'MaxIter',1e10) ;
capT        = find(Month==0) ;  
capT2   	= 11 ;
date        = linspace(1,capT2,capT2) ;
Nw          = 100;
Nwlong      = 1000 ;
b_m         = .01 ;
b_f         = .01 ;
zeta_m      = .6 ;
zeta_f      = .6 ;
mustar_m    = [mean([MeanLogAdjustedOfferedCompMen((1:capT)) ],'omitnan')] ;                  % these are the actual values
sigstar_m   = [mean([StandardDeviationLogAdjustedOfferedCompMen((1:capT)) ;StandardDeviationLogAdjustedOfferedCompWomen(1:capT)],'omitnan')];   
mustar_f    = [ mean(MeanLogAdjustedOfferedCompWomen(1:capT),'omitnan')] ;                  % these are the actual values
sigstar_f   = [mean([StandardDeviationLogAdjustedOfferedCompMen((1:capT)) ;StandardDeviationLogAdjustedOfferedCompWomen(1:capT)],'omitnan')];   

h1 = figure(1)
plot(MeanLogAdjustedOfferedCompMen,'b')
hold on
plot(MeanLogAdjustedOfferedCompWomen,'r')


h2 = figure(2)
plot(StandardDeviationLogAdjustedOfferedCompMen,'b')
hold on
plot(StandardDeviationLogAdjustedOfferedCompWomen,'r')
hold off
close all

mu0_m       = mustar_m ;
mu0_f       = mustar_f ;
lam_m       = [mean([ProbofReceivinganOfferifActivelySearching1819cohortsMen(6:capT+6); ProbofReceivinganOfferifActivelySearching1819cohortsWomen(6:capT+6)],'omitnan')] ;
lam_f       = [mean([ProbofReceivinganOfferifActivelySearching1819cohortsMen(6:capT+6); ProbofReceivinganOfferifActivelySearching1819cohortsWomen(6:capT+6)],'omitnan')] ;
gam_m       = 0 ;
gam_f       = 0 ;
iota_m      = 2 ;
iota_f      = 2 ;


